The origin of hysteresis in resistive switching in magnetite is Joule heating. 
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In many transition metal oxides the electrical resistance is observed to undergo dramatic changes 
induced by large biases. In magnetite, Fe304, below the Verwey temperature, an electric field 
driven transition to a state of lower resistance was recently found, with hysteretic current-voltage 
response. We report the results of pulsed electrical conduction measurements in epitaxial magnetite 
thin films. We show that while the high- to low-resistance transition is driven by electric field, the 
hysteresis observed in I — V curves results from Joule heating in the low resistance state. The shape 
of the hysteresis loop depends on pulse parameters, and reduces to a hysteresis-free "jump" of the 
current provided thermal relaxation is rapid compared to the time between voltage pulses. A simple 
relaxation time thermal model is proposed that captures the essentials of the hysteresis mechanism. 

PACS numbers: 71.30.-|-h,73.50.-h,72.20.Ht 



Dramatic changes in resistance induced by electric 
fields, so called resistive switching (RS), have recently 
attracted much attention due to this phenomenon's po- 
tential application in memory devices (resistive random 
access memory, ReRAM) [il, [3. RS from high- to low- 
resistance states is driven by application of high volt- 
age, and corresponding up-and-down sweeps of current- 
voltage {I-V) characteristics often show hysteresis, i.e. 
in sweeps up and down in bias voltage, the current does 
not retrace itself. Systems exhibiting hysteretic RS in- 
clude organic compounds 0| and transition-metal oxides 
such as widely-studied colossal resistance manganites 4 1 , 
perovskites (e.o. SrTiOa 0), ID cuprates Sr2Cu03 6|, 
NiO 0, TiOa Q etc. 

For some RS systems, while sweeping out a hysteresis 
loop in I-V with a switch to a low resistance state at 
high bias, the low resistance state persists down to zero 
current as voltage approaches zero. This behavior is often 
the case for RS systems where the switching is based 
on metallic filament formation at a transition point Q. 
However, for some RS systems the low-resistance state 
persists only in some voltage interval, and the system 
returns to the high-resistance state before voltage returns 
to zero. This is the case for some complex oxidesjl, 
[To| as well as for magnetite nanostructures, which were 
recently shown to exhibit RS at low temperatures [111, 

Magnetite, Fe304, is an example of strongly corre- 
lated material. In equilibrium, bulk magnetite under- 
goes a structural transition at the Verwey temperature, 
Tv ^^120 K, accompanied by three-order-of-magnitude 
change in electrical conductivity, i.e. a metal- insulator 
transition (MIT) [l^. Recently we demonstrated that 
magnetite nanoparticles and thin films, once in the in- 
sulating state below Ty, exhibit RS under a sufficiently 
large voltage bias [ll|. By examining RS systematically 
in different device geometries, the switching was demon- 



strated to be driven by the applied in-plane electric field. 
This is in contrast to previously observed transitions in 
magnetite driven b y J oule heating of the samples above 
Tv under bias [3, When the voltage is swept con- 
tinuously, the electric field-driven switching takes place 
at either polarity of voltage with well-pronounced hys- 
teresis. 

In this paper we determine the origin of hysteresis 
in magnetite structures based on epitaxial thin films. 
Through extensive voltage pulse measurements with con- 
trolled pulse parameters we unambiguously show that, 
while the high resistance to low resistance switching is 
driven by electric field, the hysteretic behavior originates 
from local Joule heating of the channel once the system 
is switched to the low-resistance state. A very simple 
thermal model agrees with the data with appropriately 
chosen parameters. The parameter values required to 
achieve quantitative consistency with the data demon- 
strate that one must consider heating and thermal trans- 
port beyond just the magnetite film itself. 

Epitaxial magnetite thin films with a 50 nm thick- 
ness were grown on (100) MgO single crystal substrates 
by oxygen-plasma-assisted molecular beam epitaxy. De- 
tails of the growth process have been reported elsewhere 
0; 113 • The films were characterized by reflection high- 
energy electron diffraction (in situ, during growing), high 
resolution X-ray diffraction measurements, Raman spec- 
troscopy and resistance measurement to prove crystalline 
quality and stoichiometry of the samples [l^, [l0| . The 
films show a jump in the temperature dependence of re- 
sistance at T ~ 108 K (Fig. [T^ inset), characteristic of 
the Verwey transition in magnetite thin films. 

Devices for two-terminal measurements were prepared 
by electron beam lithography. A channel length 190 - 
900 nm is defined by two 5-10 /im wide current leads 
(Fig. [1^ inset) which are connected to micrometer-size 
pads (300 X 300 ^m). Upon testing several different 
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contact metals (Au, Pt, Cu, Fe and Al), copper showed 
the lowest contact resistance with the magnetite film. 
Thus, the electrodes were made of 6 nm Cu / 15 nm Au 
thin layers deposited by electron beam evaporation. 

Electrical characterization of the samples was per- 
formed by a standard two-terminal method using a semi- 
conductor parameter analyzer (Hewlett Packard 4155A). 
The voltage was applied to the source lead with the drain 
grounded, and current flowing through the channel was 
monitored. Two different measurement methods were 
used: a continuous staircase sweep of the source volt- 
age up and down; and a pulsed sweep implemented with 
controlled pulse parameters. Each voltage sweep consists 
of 1000 points equally spaced in voltage. Measurements 
were performed in the temperature range 80-300 K with 
5 K steps. 

At high temperatures I-V curves have a slightly non- 
linear shape (Fig.[T^, b) typical for many transition metal 
oxides. Below a certain temperature, coincident with 
the Verwey temperature inferred from measurements of 
the zero-bias resistance, a sharp jump in current was ob- 
served as the source voltage reached a critical value, as 
described previously [ll|. At 80 K, for example, the cur- 
rent changes abruptly by a factor of 100 at a particular 
switching voltage. The state after transition ("On" state) 
shows an approximately linear I-V dependence with a 
much smaller differential resistance than that of high- 
resistance ("Off") state prior to the transition. At a 
given temperature and device geometry, the switching 
happens at a critical voltage, V^"", and as the voltage 
is swept back the system remains in the On state un- 
til it switches back to the Off state at a switch-off volt- 
age, V^ff . For a continuous staircase sweep V^^^ is 
always lower than V^, resulting in well defined hystere- 
sis (Fig.[T^). The transition is symmetrical along V-axis, 
with identical transitions occurring at positive and neg- 
ative voltage sweeps. 

Switching phenomena are observed at temperatures 
right below Verwey temperature ('^108 K) and, thus, 
magnetite is comparatively insulating in the undisturbed 
state (see zero-bias RvsT dependence in Fig.[T^ inset). 
Previously we demonstrated that the observed transition 
is driven by electric field, not by thermal heating of the 
magnetite sample above Verwey temperature 11[ . In the 
high resistance state (1^1 < V^"") heating of the sample 
is comparatively negligible, whereas in the On state dis- 
sipated power significantly increases, making Joule heat- 
ing of the sample much more likely. Below we use pulsed 
measurements to demonstrate that this On state heating 
is the origin of the apparent hysteretic behavior in our 
system. 

We performed pulsed sweep experiments, with con- 
trolled duration of the voltage pulse (pulse width), the 
time between consecutive pulses (pulse period) and the 
voltage value at rest between the pulses (base voltage) 
(Fig. [T]3 inset) . In these pulsed measurements the base 
voltage may be kept Vbase = V, eliminating Joule heat- 
ing and allowing cooling of the channel between pulses. 
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FIG. 1: Typical I-V curves at different temperatures in con- 
tinuous staircase (a) and pulse (b) sweep modes. Arrows indi- 
cate the direction of voltage sweeps. Above a certain temper- 
ature (in this case ~ 105 K) sharp jumps in current are not 
observable. Top (a) inset shows a scanning electron micro- 
graph of a two-terminal device on the surface of a magnetite 
film; bottom (a) inset is a temperature dependence of resis- 
tance demonstrating Verwey transition at Tv ~ 108 K. The 
(b) inset schematically represents the pulse parameters for the 
pulsed sweeps. Successive applied voltages differ by 1-15 mV 
depending on the total sweep range. 

The dependence of system response on pulse width and 
pulse period can indicate heating dynamics of the sample 
(especially in the On state). 

First, we investigated the system response to pulse 
sweeps with different widths: we varied the pulse duration 
while keeping fixed the time system stays at base voltage 
(pulse period - pulse width = const, {Vbase = V), see 
fig- Hb inset). Fig. [5^ shows resultant I-V curves (only 
positive voltages for better visualization) at 80 K. The 
shape of the hysteresis clearly changes as pulse width 
increases: V^^^ moves farther away from V^, which 
remains independent of the pulse width (Fig. [5^). Note 
that the current is consistently higher for longer pulses. 

In a second set of pulse experiments we made volt- 
age sweeps with various pulse periods and constant pulse 
width. The shortest possible pulse width (0.5 ms in our 
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FIG. 2: Dependence of hysteresis (J-V curve) shape at 80K 
in pulse sweep mode on (a) pulse width while the time sample 
rests at Vtase is kept constant (4 ms) and (b) pulse period, 
while pulse width is kept the same (1 ms). Base voltage, Vtase, 
is always set to V. 

setup) was used to minimize heating of the sample. The 
increase of pulse period has an opposite effect on V^^^ 
position than the increase of pulse width: V.'if^ moves 
closer to as the pulse period increases (Fig.[2lD), i.e. 
as the system rests at zero bias voltage between pulses for 
longer times. Note, again, that position remains in- 
dependent of pulse period. At long enough period times 
(> 100 ms) V^" and ^ are the same value. Fig. [T)d 
shows I-V curves at different temperatures obtained by 
pulse sweeps with a 100 ms pulse period. Instead of hys- 
tereses there are sharp jumps in On state and back to Off 
state at almost the same voltage: V^J^ ~ ■ 

While hysteretic RS (with differentiated V^j^ and 
V2J^ values) was previously observed in many systems, 
here we demonstrate that in magnetite there is a single 
intrinsic switching value, Vsw = = ^sv/^ ^ ^^'^ 

intrinsic hysteresis. Thus, the switching voltage is a char- 
acteristic parameter of the system: at every temperature 
and device geometry, there is a certain switching volt- 
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FIG. 3: Temperature dependence of switching voltage (open 
squares) and its exponential fit (solid line). The inset shows 
I-V curves in the On state in pulse sweep mode with the 
same pulse width (0.5 ms) but difi^erent pulse periods, showing 
that rapid pulse repitition leads to more severe heating, as 
expected. 

age required to drive the transition. This K,to decreases 
as temperature approaches Ty, following a nearly expo- 
nential dependence (Fig. [3]) . The separation of and 
V2J^ values, shifting of V^J^ position in various pulse 
period and width sweeps, and the observation of hystere- 
sis are easily explained by heating of the sample in the 
On state. 

We first give a qualitative picture of how On state heat- 
ing leads to hysteresis, and then provide a quantitative 
description with a simple model of channel temperature 
during a pulse sweep. The sample begins at an initial 
temperature Tset \ for this temperature and a given chan- 
nel length, there is a corresponding switching voltage, 
yswiTset) (Fig. [3]). Consider performing a sweep with 
pulse parameters selected to give hysteresis. First, as we 
start sweeping voltage from V up, the system remains 
approximately at Ts^t prior to the transition point due to 
the comparatively small amount of Joule heating. After 
application of a pulse with V > Vsw (Tset), the system 
is driven into the On state. With the resulting increased 
conductance (and, thus, current) the Joule heating of the 
channel is much larger, elevating the local effective tem- 
perature of the channel to Tcur > Tset ■ Non-stcady-state 
warming of the sample continues further in the On state. 
This is clear from the non-linear increase of the current 
and the fact that the I-V trace in the On state does 
not retrace itself (Fig. [3] inset). Note that sufficiently in- 
creasing the pulse period does allow I-V curves to retrace 
themselves on up and down sweeps (Fig. [3] inset). As the 
pulse voltage is swept back below Vsw {Tset), the channel 
remains in the On state because the channel is at the 
elevated effective temperature, Tcur, and the switching 
voltage is therefore lower: VswiTcur) < VswiTset)- 

If the channel cools down to Tset between two consec- 
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utive pulses, the system remain in the OS state for any 
voltage applied that is lower than VgwiTsgt) and jumps 
are observed instead of hysteresis (Fig. [Ijj). If the pulse 
period is short compared to the thermal relaxation time, 
then the channel remains at Tcur > Tset, and we ob- 
serve the separation of = Vsw{Tset) and V^^^ = 
VswiTcur), and resulting hysteresis. Upon increasing the 
pulse period (system stays longer at Vbase = V), the 
channel has more time to cool after application of previ- 
ous pulses, Tcur is closer to Tget and V^^^ moves closer to 
in excellent agreement with the experimental data 
(Fig. [^b)- In experiments on increasing pulse widths, it 
is clear that longer pulses heat the channel more than 
shorter pulses, thus Tcur is higher in the former case, 
corresponding to lower VswiTcur)- Thus, experimentally 
observed moving V^^^ away from V^J^ as the pulse width 
increases (Fig. [2^) is completely consistent with this heat- 
ing picture. 

As qualitatively described above, V^^^ position is re- 
lated to the current local temperature of the channel. To 
quantitatively describe the movement of ^ position 
in sweeps with different pulse parameters, we should con- 
sider a model to estimate the effective channel tempera- 
ture as a function of applied pulsed voltage (and time). 
Detailed thermal modeling is very challenging because of 
the highly local character of the heating and the diffi- 
culty of capturing all the relevant heat transfer processes 
in these nanostructures. Instead we consider a simple 
model where the return to equilibrium state (relaxation) 
is described by a relaxation time, r (20| . and the temper- 
ature change is: 
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where T is a current temperature of the channel. The 
heat flux into the sample is obviously dQin/dt = V x I. 
We further assume the heat dissipation can be simply 
described as 



dQo 



dt 



T-n. 



(2) 



where Cy is a heat capacity of the sample (in Joule/K). 

Note that this model assumes Cy and r to be inde- 
pendent of temperature, which is clearly an idealiza- 
tion. The temperature variation with respect to time, 
dT/dt, considering total heat flux through the system, 
i.e. §-AQi7i + Qout), is given by: 



Qout) ^Vxl-Cy 



T-T,, 



(3) 



By solving Eq. ([3]) we derive an expression for a temper- 
ature, T(t + dt), given the current temperature, T(t), in 
a moment of time, dt: 



T{t + dt) = T„ 



V X / X r , , dt . . 
X (1 -exp( ))- 
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FIG. 4: Calculated temperature in relaxation time model (see 
Eq. (|4])) with respect to source voltage (a) and correspond- 
ing time (b) for a 3 ms-wide / 7 ms-period voltage pulse 
sweep. Black points reflect T values before the transition 
point, V < Vsw After the transition point, blue triangles 
show T immediately after pulse application and red circles 
show T after relaxation between two sequential pulses. The 
hatched area in (b) indicates the time interval depicted in 
detail in (c): Applied pulsed voltage as a function of time 
(top) with corresponding calculated temperature of the chan- 
nel (bottom). Shaded areas indicate time intervals when sys- 
tem is is not under bias {i.e. source voltage is zero) and only 
relaxation to Tset is taking place (see Eq. ((Sjl). 



iT{t)-Tset) xexp(--) 

T 



(4) 



When V ^ 0, both heating and cooling (by relaxation 
to equilibrium Tget) processes occur described by second 
and third terms of the right side of Eq. ^ , respectively. 
Between pulses (Vf,ase=0 V) only heat dissipation takes 
place and Eq. ^ reduces to the expected exponential 
decay of the temperature with the relaxation time, r, as 
a scaling factor: 

Tit + dt) = T,ct + iT{t) - T,ct) X exp (--) (5) 

T 

There are two fitting parameters in this model that gov- 
ern the evolution of the current temperature: heat ca- 
pacity, Cy, and relaxation time, r. 

An example of T variation calculated by Eq. ([4]) is 
shown in Fig. 3] with real experimental data for sourccd 
voltage, V , and measured current, /, (3 ms-width, 7 ms- 
pcriod pulse sweep) and Cy and r values as will be 
discussed below. Figs. |4^, b show the calculated tem- 
perature of the channel as a function of source voltage 
and corresponding time, respectively, over a large volt- 
age range around transition point. Black points represent 
T values before the transition to the On state and they 
do not deviate much from Tset (80 K in this case). As 
the system is switched to the ON state the temperature 
profile changes drastically: As pulse voltage is sourced. 
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T increases and blue triangles show the calculated T val- 
ues right after pulse application ("high-T" state). After 
relaxation for a time between pulses, the channel cools 
down to temperatures represented by red circles ( "low-T" 
state) which are only slightly higher than original T^et- 
For better visualization, Fig. zooms into the small re- 
gion indicated by the hatched area in Fig. at a turning 
point of a sweep (from sweep up to sweep down in volt- 
age). It clearly demonstrates how the temperature first 
rapidly rises as the pulse is applied and then decreases as 
the system rests at Vbase = OV^ (shaded regions), relax- 
ing back toward Tset = 80 K. Both "high-T" and "low-T" 
temperatures follow the applied voltage: T increases as 
V is swept up and starts decreasing as sweep is reversed 
and V decreases (follow the red circles and blue rectan- 
gulars in Fig. Hj: bottom). 

Each device has its defined single- valued Vsw (T) func- 
tion (an example is shown in Fig. [3|) with Vsiu being a 
characteristic parameter at a certain temperature. Dur- 
ing a voltage sweep down, the temperature at the mo- 
ment the pulse is applied decreases (follow the red circles 
in Fig. [3):) while Vsw corresponding to this temperature 
increases. Note that when the pulse voltage is applied 
the relevant temperature is that of the "low-T" state (red 
circles in Fig.[3t bottom), thermally relaxed after appli- 
cation of the preceding pulse. Switching from On to Off 
states happens as soon as applied pulse V value is lower 
than Vsw for the current temperature of the channel at 
the time the pulse is applied. 

We find that the proposed model can match experi- 
mental V^^^ positions with appropriately chosen Cy and 
T values. For each device at a given Tset we can find a 
Cy and T pair which accurately describes the shifting of 
V^f^ position for sweeps with different pulse widths and 
the same time between pulses (pulse period - pulse width 
= const). An example of temperature dependence on ap- 
plied pulsed voltage is shown in Fig. [5] corresponding to 
experimental I-V curves in Fig. [2^. The adjusted opti- 
mum values of C„ and r are 10~^ Joule/K and 1.5x10^'^ 
seconds, respectively. 

To get a sense of what these parameter values imply, 
we note that the estimated heat capacity value based on 
specific heat capacity and density of magnetite {C„ ~ 40 
Joule/(molxK) at 80 K [HI, p = 5.18 g/cm^) and geo- 
metrical parameters of the channel is ~ lO^^""* Joule/K. 
This is much lower than Cy value required by the model 
to match the observed trends, strongly implying that 
heating is not confined only to the channel, but involves 
a much larger volume. This also implies that the MgO 
substrate and Au electrodes are relevant to the thermal 
relaxation process. 

For sweeps with different pulse periods and the same 
pulse width, VsJ,^ positions can not be satisfactory mod- 
eled with a single Cy and r pair. For a fixed C„, 
sweeps with longer pulse periods demand larger r val- 
ues to match the experimentally observed V°Jj^ positions. 
Longer period times imply larger temperature variation 
over the time, relaxation process takes place (sec Eq. ([S])). 
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FIG. 5; The dependence of calculated temperature on source 
voltage in the relaxation time model for sweeps with different 
pulse widths and the same time between pulses (4 ms), Tset 
= 80 K. "High-T" (immediately after pulse application) and 
"low-T" (after relaxation between pulses) state temperatures 
are plotted in the same color, with the latter apparent at the 
bottom of the figure. Arrows indicate switching back position 
for each sweep which happens when applied voltage becomes 
lower than Vsw for a current temperature of the channel. 

This is almost certainly due to the failure of the ideal- 
ized assumption of Cy and t being independent of tem- 
perature. A much more sophisticated thermal model (in- 
corporating the spatial temperature distribution and the 
temperature dependent thermal properties of the mag- 
netite, the substrate, and the electrodes) is likely neces- 
sary for a complete quantitative picture. 

Since heating effects in the On state are responsible for 
hysteresis in this system, it is worthwhile to review the 
arguments that the high resistance to low resistance tran- 
sition itself is NOT driven by heating the sample above 
Verwey temperature, but indeed is driven by electric field 
[ll[. First, a strong argument against heating effect is the 
temperature dependence of the power, Psw, dissipated at 
the transition point. It was shown in Ref. [Tl| for two- 
terminal devices that P^^ = V x I at V^ decreases 
as temperature decreases. This behavior is incompati- 
ble with the assumption that raising of the channel tem- 
perature above Verwey temperature (here Ty ^108 K) 
is responsible of the resistance switching. Assuming a 
heating effect, the lower the set temperature the more 
power should be provided to warm up a channel above 
Ty. the opposite of the trend observed in the experiment. 
In this work we confirm the Psw vs T dependence in the 
two-terminal geometry (Fig. [6]). Additionally, we per- 
formed four-terminal measurements (details to be given 
separately [22j) and calculated the corresponding local 
switching power dissipated in the channel, Pfw = x / 
at a transition point, where AV is a voltage difference 
between two voltage probes within the channel. P^^ is 
much lower than P^!^ because of the significant contact 
voltage drop at the electrode/Fe304 interface. The power 
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FIG. 6: Temperature dependence of the power at a transition 
point, , in two-terminal (PI"^ = V x I) and four-terminal 
(Paw ~ X experiments. Inset shows the dependence 
of the two-terminal switching voltage, Vs™, on the channel 
length at two different temperatures. Solid lines represent 
a linear fit; the slope of each line reflects the electric fleld 
strength required to drive the transition. 

dissipated directly in the channel, Pj'J, shows the same 
temperature dependence as P'^'^ (Fig. (H) , arguing further 
against simple heating driving the transition. 

A second argument for that the transition is electri- 
cally driven is the dependence of Ysw on the length of the 
channel. At each temperature V^w scales linearly with the 
channel length, L (Fig.[n]inset). This implies that at each 
temperature there is a certain electric field necessary to 
drive a transition. Extracted from the slope of Yaw vs. 
L, the critical electric field value is about 7x10"' V/cm 
at 105 K and it increases slightly up to ~10^ V/cm down 
to 80 K (Fig. [6] inset). These values are much lower than 
the catastrophic dielectric breakdown field (~10^ V/cm) 
[23! for most insulators. 

Third, pulse experiments that widely span heating con- 
ditions show no shift in K,^ position upon varying cither 
pulse width or pulse period. Longer pulses deliver more 



energy to the channel (= / x y x where t is the dura- 
tion (width) of the pulse) . Assuming a thermal switching 
mechanism, the switching would be expected at lower 
voltages for longer pulses, which is not the case in the 
experiments. Different pulse widths and periods have 
tremendous effects (both in V^J^ positions and current 
values) on the behavior of the system in the On state 
(Fig. (2) when heating is clearly an issue. In the Off state 
only minor changes in 1-V curves were detectable for 
sweeps with all available widths and periods. Finally, 
sourcing a constant voltage for minutes right below Vgw 
value at a certain temperature does not induce or alter 
the transition. 

In conclusion, resistive switching in magnetite thin 
films from a high resistance state to a low resistance state 
is driven by electric field, while the hysteresis observed in 
such switching is a result of thermal effects. We find that 
a simple relaxation time model of the thermal processes 
can describe the shifting of V°^^ positions for sweeps 
with different pulse parameters. Model parameters indi- 
cate that the substrate and electrodes are important in 
determining the thermal dynamics of the system. The 
electric field-driven switching observed in magnetite is 
different than several mechanisms suggested for RS in 
perovskites, such as dynamics of oxygen vacancies and 
interfacial effects. First, the lack of intrinsic hysteresis 
differentiates magnetite RS from that in the perovskites. 
Second, magnetite RS can be induced only below the Ver- 
wey temperature, Ty, whereas RS in perovskites is ob- 
served in a wide range of temperatures [al . Below Ty 
magnetite is in a correlated ordered state [2^. [25t. Thus, 
RS in Fe304 is a strong candidate for the theoretically 
predicted breakdown of charge- ordered states by electric 
field p6| . The mechanism of this nonequilibrium transi- 
tion demands further investigation. 
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